The zero-energy state in graphene in a high magnetic field 
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The fate of the charge-neutral Dirac point in graphene in a high magnetic field H has been 
investigated at low temperatures (T ~ 0.3 K). In samples with small gate- voltage offset Vq, the 
resistance _Ro at the Dirac point diverges steeply with H, signalling a crossover to an insulating 
state in high field. The approach to the insulating state is highly unusual. Despite the steep 
divergence in Rq, the profile of Rq vs. T in fixed H saturates to a T-independent value below 2 K, 
consistent with gapless charge-carrying excitations. 

PACS numbers: 73.63.-b,73.21.-b,73.43.-f 



The discovery of the quantum Hall effect (QHE) in 
monolayer graphene crystals provides a new system for 
investigating relativistic Dirac-like excitations in solids 
0, 0, U0, In a magnetic field H, the system forms 

Landau Levels (indexed by n) that are 4-fold degenerate. 
The Hall conductivity a xy is accurately quantized as the 
chemical potential /x is changed from the hole part to elec- 
tron part of the Dirac spectrum. Considerable attention 
has focussed on the n = Landau Level (LL), especially 
on the nature of the electronic state at the charge-neutral 
point (/i = 0) in an intense magnetic field H. Several 
groups 0, S & 01, ED, El have predicted that a large 
field stabilizes the quantum Hall ferromagnetic (QHF) 
state, in which the pseudospins describing the valley de- 
gree of freedom become ferromagnetically ordered. In a 
second group of theories [HI, Il4j , interaction causes an 
excitonic gap to open at the Dirac point. Experiments 
are actively addressing these issues [HI, El [TtJ • Jiang et 
al. [TU have inferred that the sublevel gaps at v = and 
±1 arise from lifting of the spin and sub-lattice degen- 
eracies, respectively, and inferred a many-body origin for 
the states. We have found that, in samples with small Vq 
(the gate voltage needed to align \i with the Dirac point), 
the value Rq of the resistance R xx at the Dirac point di- 
verges steeply with H, i.e. a large H drives the Dirac 
point insulating. Despite the strong H dependence, Rq 
saturates to a T-independent value below 2 K, provid- 
ing evidence for charged, gapless excitations. In samples 
with large Vq, this divergence in Rq is shifted to higher 
fields. 

Following Refs. 0, H, 0|, we peeled single-layer 
graphene crystals (3-10 fxm in length) from Kish graphite 
on a Si-Si02 wafer. Au/Cr contacts were deposited us- 
ing e-beam lithography (Fig. [Th., inset). We have found 
that the high-field behavior of Rq is strongly correlated 
with Vq (Table H]). All samples (except K22) have fj, ly- 
ing in the electron band (positive Vq). Samples in which 
| Vol < 1 V (K7 and K22) display a very large R a (U) 
(resistance measured at 14 T and 0.3 K), which arises 
from the strong divergence mentioned. By contrast, in 
samples with large |V |, i2o(14) < 7 kfi. 

Figure [T^l shows the variation of R xx in K7 plotted vs. 
the shifted gate voltage V' = V g — Vq with H held at 8, 11 
and 14 T (at T = 0.3 K). The striking feature here is that 



the peak corresponding to the n — LL increases to >100 
kO at 14 T, whereas the peaks corresponding to n = ±1 
remain below -7 kfl. As in Refs. [J 0, H II II, the 
Hall conductivity a xy (Panel b) displays plateaus given 
by El 
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where n indexes the 4-fold degenerate LL and v indexes 
individual sublevels. In K7, the 'zero' plateau a xy ~ at 
Vg = is already visible at H = 8 T. 

Narrowing our focus to the n = LL, we examine 
R xx in the n = LL as a function of T, with H fixed 
at 14 T (Fig. [2}i). We see that, from 40 to 0.3 K, Rq 
rises steeply from 4 kil to 190 kfl. The curve of the 
conductivity <r xx plotted vs. V' reveals a 2-peak struc- 
ture that implies splitting of the 4-fold degeneracy by a 
gap A (Fig. [2p). At 100 K, the two peaks are already 
resolved. With decreasing T, the minimum at Vg = 
initially deepens rapidly, but saturates below ~2 K. The 
Hall conductivity a xy at 0.3 K (thin curve) displays a 
well-defined plateau on which a xy ~ 0. Because the next 
plateau is at a xy = 2(e 2 /h), we infer that each of the 
peaks in a xx is comprised of 2 unresolved sublevels. The 
opening of the gap causes a xx (at V' g — 0) to fall rapidly 
with decreasing T, until saturation occurs below 2 K. 

The behavior of Rq described here is qualitatively dif- 
ferent from that in, for e.g., Ref. [l5j . To understand the 
difference, we have examined several samples (Table H]) . 
As mentioned, the offset gate Vq is a crucial parameter. 
Figure compares the curves of R xx (n = LL) in the 



Sample 


Vb(V) 


ifc(14)(kn) 


Ml/T) 


K5 


3 


80 


0.3 


K7 


1 


190 


1.3 


K8 


12 


15 


0.6 


K18 


20 


7.5 


0.9 


K22 


-0.6 


>280 


2.5 


K29 


22.5 


7 


0.2 



TABLE I: Sample parameters. Vo is the gate voltage needed 
to bring /x to 0. 7?o(14) is Ro measured at H = 14 T and T 
= 0.3 K. fie is the electron mobility at H — 0. 
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FIG. 1: (color online) The resistance R xx (Panel a) and Hall 
conductivity a xy (b) in Sample K7 versus (shifted) gate volt- 
age Vg = V B - Vo at 0.3 K with H fixed at 8, 11 and 14 T. 
Peaks of R xx at finite Vg correspond to the filling of the n = 1 
and n — 2 LLs. At Vg = 0, the peak in R xx grows to 190 
kfi at 14 T. The inset shows in false color a graphene crystal 
(dark red) with Au leads deposited (yellow regions) . The bar 
indicates 5 /xm. Panel b shows the quantization of a xy at the 



values (Ae 2 /h)(n+\] 
around V g ' = 0. 
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samples K5, K7, K8 and K29, all measured at 14 T at 
T = 0.3 K. For each sample, we have plotted R xx vs. 
the unshifted gate voltage V g , so its peak automatically 
locates Vq . It is clear that K7 (Vo — 1 V) has the highest 
peak, followed by K5 (V = 3 V), whereas K8 (V = 12 
V) and K29 (22.5 V) have peaks that are severely sup- 
pressed. Further insight into this pattern of suppression 
is given below when we examine the H dependence of Rq 
at low T. 

Sample self-heating may also obscure the divergence. 
We find that, below 1 K, self-heating becomes serious 
when the dissipation exceeds ^2 pW. The measurements 
of R xx vs. V g were repeated at 3 currents (I— 0.6, 2 
and 15 nA) at T = 0.3 K. The results at I = 0.6 and 
2 nA are virtually identical. However, the curve at 15 
nA is 30% smaller near V' — 0, consistent with heating. 
Hence, we have kept I at 2 nA to eliminate self-heating as 
a problem. Heating at the contacts is negligible because 
of the small contact resistances (<~1 kfi) relative to Rq. 

Hereafter, we focus on Rq, or equivalently, the Dirac- 
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FIG. 2: (color online) The resistance R xx , conductivity o xx 
and the Hall conductivity tj xy in K7 vs. the shifted gate volt- 
age Vg, with H fixed at 14 T. As T decreases to 0.3 K, the 
zero-energy peak in R xx (Panel a) rises steeply to 190 kfi. 
Panel (b) shows that, as T decreases, double peaks in o xx 
are clearly resolved. Between the peaks, o xx falls rapidly but 
saturates below 2 K. The Hall conductivity at 0.3 K (thin 
curve) displays a clear plateau (\(J xy \ < 0.02e 2 //i) in the in- 
terval — IV < Vg < IV. Panel (c) compares R xx (of n — 
LL) vs. unshifted gate V a in the samples K5, K7, K8 and K29 
at 0.3 K. In each sample, R xx peaks at Vo. As Vo increases, 
i?o(14) rapidly decreases. 



point conductivity a xx = L/wRq (L and w are the length 
and width). Curves of the conductivity versus log 10 T are 
shown in Fig. [3^ at selected fields. In low fields (H < 
9 T), the T dependence of <j xx is quite mild. As H is 
increased to 14 T, the opening of the gap A (between 
the n=0 sublevels) causes the conductance to decrease 
sharply below 40 K. However, instead of falling to 0, a xx 
saturates below 2 K to a T-independent residual value 
a res , as anticipated in the discussion of Fig. ^jp. The 
existence of this residual a res , which is highly sensitive 
to H , is one of our key findings. 

The field dependence of a res is best viewed as a diver- 
gent Rq. Figure [3Jd shows the rising profile of Rq vs. H 
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FIG. 3: (color online) The T dependence of a^ x in K7 (= 
L/wRo) and the H dependence of Ro at low temperature. 
Panel a shows curves of a® x vs. log 10 T with H fixed at 6-14 
T. For H > 8 T, the gap A causes a® x to decrease markedly 
until saturation at the residual value a re3 occurs below 2 K. 
Panel (b) displays the steep increase in Ro vs. H in K7 at 
selected T. At 0.3 K, Ro appears to diverge at a field near 
18 T (see Fig. |4p). Panel (c) compares the Ro{H) profiles in 
Samples K7, K18 and K22. In Sample K18 (Vb = 20 V), the 
divergence in Ro becomes apparent only above 14 T, whereas 
in K22 (Vo = -0.6 V) Ro starts to diverge at fields lower than in 
K7. In K7, we have plotted Ro values measured by sweeping 
V g at fixed H (solid symbols) with Ro measured by sweeping 
H with Vg fixed at (solid curve), to show consistency. 



in sample K7 at selected temperatures. The divergent 
form at the lowest T (0.3 K) strongly suggests that the 
system is rapidly approaching a field-induced crossover 
(or transition) to an insulating state. 

In light of the importance of Vb, it is instructive to see 
how the profile of Ro vs. H varies between samples. Fig- 
ure [3J; compares the results in K7, K18 and K22 at T = 
0.3 K. In K18, where Vb (20 V) is quite large, the diver- 
gence in Rq(H) becomes noticeable only in fields above 
14 T. Conversely, in K22 for which V (-0.6 V) is slightly 



FIG. 4: (color online) (Panel a) The contour plot of Ro(T, H) 
(K7) in the T-H plane (vertical bar shows values of Ro). The 
contour lines emphasize the unusual approach to the insulat- 
ing state. At low T, Ro is unchanged on a horizontal path (H 
held fixed) , but it rises rapidly on a vertical path (increasing 
H at fixed T). Panel (b) displays logJ?o vs. l/%/l — h in 
K7 at T = 0.3 K, with h = H/Ho, where H = 18 T. The 
linear segment at large Ro shows that the divergence is consis- 
tent with Ro(h) ~ exp[2fe/Vl - h] with b ~0.7. In Panel (c), 
the plot of {din Ro / dH)~ 2 / z vs. H shows a high-field linear 
segment that extrapolates to zero at ~ Ho (18 T). 



smaller than in K7, Rq diverges at field scales smaller 
than in K7. From the trend, it is clear that the diver- 
gence in i? is shifted to ever higher fields as Vb increases. 
Referring back to Fig. [^fc, we now see that the strong 
suppression of i?o(14) in samples with large Vb simply 
reflects the shift of the divergence to larger H. These re- 
sults underscore the importance of choosing samples with 
| Vb| < 1 V for investigating the intrinsic properties of the 
Dirac-point. In the case of K7, we have also checked that 
the 2 procedures for measuring Ro (varying V g at fixed 
H or vice versa) give consistent values. 

The variation of Ro (T, H ) in K7 is conveniently rep- 
resented in a contour plot in the T-H plane (Fig. 
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Below ~2 K, the contour lines are horizontal, which im- 
plies that i?o is unchanged if the sample is cooled in fixed 
H. This provides evidence that a res involves gapless ex- 
citations. However, if T is fixed, Rq rises steeply with 
H , implying proximity to the insulating state (deep-red 
region). When a system approaches the insulating state, 
its resistivity generally diverges as T — ► 0, as a result 
of either strong localization (variable-range hopping) or 
the opening of a mobility gap (weak localization is not 
relevant here because of the intense H). In both cases, 
decreasing T reduces the conductance because the itiner- 
ant states are severely depopulated. Hence, the pattern 
in Fig. [4^ is most unusual. The gaplessness of a res sug- 
gests that, below 2 K, these excitations are protected 
from the effects of changing T. Paradoxically, they are 
not protected from an increasing H, which reduces the 
current carried at an ex pon ential rate. 

In the theory in Refs. [l5lll^|. the current at the Dirac 
point is carried by a pair of edge states. The change in 
R xx vs. H is interpreted as an increased scattering rate 
as the edge states are pushed closer to the edge. 

In samples with small Vq, however, the steep increase 
in i? to ^200 kQ (Fig. (3t>) suggests a different regime 
in which Coulomb exchange may be the dominant energy 
scale. From Figs. [3Jd and [4^, we infer that Ro appears 
to be diverging towards an insulating state in high fields. 
These considerations lead us to quantify the divergence 
in Rq. 

We find that Ro fits very well to the form Rq ~ 
£(/i) 2 , where the correlation length £ has the Kosterlitz- 
Thouless (KT) form 

i K T ~ exp(6/VT _ ~h), (h = H/H ), (2) 

with H replacing T. Plotting lni?o vs. vl — h, we find 
that the high-field portion becomes linear (Fig. HJa) when 
Hq is adjusted to be 17-18 T. From the slope, we find that 
the parameter b has the value ~0.7, in agreement with 
simulations of the KT transition. For self consistency, we 
may also let the data determine Hq by linear extrapola- 



tion. By Eq. H we have dlnR /dH ~ (H - H)- 3 / 2 . 
Hence a plot of (din R /dH)~ 2 / 3 vs. H should cross the 
-ff-axis at Hq. Indeed, this quantity, plotted in Fig. [4J;, 
becomes linear at large H and extrapolates to zero at 
~18 T in agreement with (b). 

Although consistency with Eq. [2] alone does not prove 
that a KT transition occurs at Hq, the fit does reveal the 
striking exponential character of the divergence in Rq. 
We adopt as a working hypothesis that this reflects the 
approach to a KT transition. The ordered state at large 
H is destroyed by the spontaneous appearance of defects 
which increase exponentially in density at H . The steep 
fall of Rq below Hq may reflect the current carried by 
these defects. 

In both the QHF state 0,1,1, El EH and the exciton- 
gap state [HI, EH , the 4- fold degeneracy in the n — LL 
is completely lifted in large H to produce an insulator. In 
the QHF, the starting symmetry is ST/ (4) if we treat the 
spin and pseudospin degrees on equal footing Q. The 
reduction of the St/ (4) to lower symmetries by the ef- 
fects of Zeeman energy, disorder or lattice discretization 
is discussed in Refs. @, E3, EH, E2- According to Ref. 
EH random disorder may force the pseudospins into the 
plane. If the ordered state indeed has U(l) symmetry, 
the XY order is susceptible to a KT transition 0, [Til ] - 
It would be very interesting to relate the fit of Rq to 
Eq. [3 with charged, topological excitations envisioned 
in the KT transition. Measurements are in progress at 
higher H to provide further evidence for a transition to 
the insulating state and to clarify its nature. 
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